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ABSTRACT 

Context. Recent analysis of stellar populations in barred galaxies have focused on the spatial distribution of stellar population ages and 
metallicities. However, barred galaxies are complex objects where dynamical instabilities play a leading role in shaping any spatial 
distribution. 

Aims. The age distribution of stellar populations should thus be analyzed from the two points of view of stellar population evolution 
and dynamical secular evolution. 

Methods. Chemodynamical simulations of single barred galaxies with simple but realistic star formation and feedback recipes are 
used to produce face-on mass-weighted maps of stellar population ages. Luminosity-weighted maps in V-band are also displayed after 
calibrating the simulation with mass-to-light ratios provided by a synthesis population model. 

Results. It is shown that inside a stellar bar two persistent diametrically opposed regions display a mean age lower than the surrounding 
average. These two low-age regions are due to the accumulation of young stellar populations trapped on elliptical-like orbits along 
the bar, near the ultra-harmonic resonance. Age gradients along the bar major-axis are comparable to recent observations. Another 
low-age region is the outer ring located near the Outer Lindblad Resonance, but the presence and strength of this structure is very 
time-dependent. 
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1. Introduction 

Thanks to the improvement in stellar libraries, population syn- 
thesis models, and inversion algorithms, dating of stellar popu- 
lations in galaxies and determination of star formation history 
has become a very active field. Indeed, the integrated spectra of 
galaxies provides a huge number of information on the current 
stellar populations and the history of their formation. Moreover, 
they simultaneously provide an insight into their current dynam- 
ical state that results both from the history of the mass assem- 
bling and their secular evolution. However, when the spectral in- 
formation is combined with the spatial one, as for integral field 
spectrographs or even long-slit spectroscopy of nearby galaxies, 
one hits a snag because the spatial variation of the star forma- 
tion history cannot be understand without any knowledge on the 
stars and gas dynamics. It is thus difficult or even impossible 
to unravel the two histories, the gravitational dynamics from the 
star formation, since the physical mechanisms at play during any 
galaxy life are so much intricately mixed. A direct approach, us- 
ing N-body and hydrodynamical simulations coupled with small 
scale recipes of star formation and feedback, is thus useful to 
disentangle the two histories even at the cost of very simplified 
physical assumptions. 

In the zoo of galaxies, barred galaxies deserve a special at- 
tention since they represent more than 75% of all disc galaxies 
in the local Universe. Most of them are still forming stars so 
that in the same object we can find several generations of stellar 
populations, of various ages and metallicities. The interpretation 
of age and/or metallicity gradients in barred galaxies are thus 



not straightforward. Using numerical simulations, Friedli et al. 
dl994|) found that a bar is able to flatten abundance gradients in 
the disc, in agreement with most observations of O/H gradients. 
They did not investigate the effects of the bar on the age distri- 
bution. 

The main goal of this Letter is thus to provide quali- 
tatively the age distribution of stellar populations in typical 
barred galaxies. Chemodynamical simulations, calibrated using 
Michel-Dansac & Wozniak (2004 ) technique, have been used to 
produce luminosity- and mass-weighted maps of stellar popula- 
tion ages. I discuss the origin and persistence of some features 
in the framework of the orbital structure of stellar bars. 

2. A simple and generic chemodynamical model 

An initial stellar population is set up to reproduce a typical disc 
galaxy. Positions and velocities for 2.5 10 6 particles are drawn 
from a superposition of two axisymmetrical Miyamoto & Nagai 
d 1975b discs of mass respectively 10 10 and 1O U M , of scale 
lengths respectively 1 and 3.5 kpc and common scale-height of 
0.5 kpc. Initial velocity dispersions are computed solving nu- 
merically the Jeans equations. The initial disc radius is 30 kpc. 
The gaseous component is represented by 50 000 particles for 
a total mass of 1.1 10 10 M0 distributed in a 6 kpc scale-length 
Miyamoto-Nagai disc. 

Evolution is computed with a particle-mesh N-body code 
which includes stars, gas and recipes to simulate star formation. 
The broad outlines of the code are the following: the gravita- 
tional forces are computed with a particle-mesh method using a 
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3D log-polar grid with (N R , N^, N z ) = (60, 64, 312) active cells. 
The smallest radial cell in the central region is 36 pc large and 
the vertical sampling is 50 pc. The extent of the mesh is 100 kpc 
in radius and +7.8 kpc in height. The hydrodynamical equations 
are solved using the SPH technique. Since I used a polar grid, 
the pre-computation of self-forces have been improved by sub- 
dividing each cell in (n r , n$, n z ) = (32, 6, 6) sub-cells. Self-forces 
are then linearly interpolated before being subtracted from grav- 
itational forces. The spatial and forces resolutions are thus much 
higher than in our previous studies based on the same code (e.g. 
Michel-Dansac & Wozniak, [2004J |2006J . 

The star formation process is based on Toomre's criterion 
for the radial instability of gaseous discs (cf. Michel-Dansac & 
Wozniak 2004 for more details). When star formation is active, 
the radiative cooling of the gas is computed assuming a solar 
metallicity. 

This simulation has been chosen as being illustrative of a 
moderately barred galaxy. 

At the end of the simulation (f*3 Gyr), the total number of 
particles is roughly 3.2 10 6 for the stellar component and 30 000 
for the gaseous one. 46% of the gas have been transformed 
into stellar particles, mainly in the central 10 kpc. However, in 
this relatively wide region, the new population (called 'young' 
hereafter) accounts for roughly 5% of the total mass (gas and 
stars) while in the central kpc, the young population accounts 
for roughly 17% of the total mass. 

This paper being focused on the expected age distribution, 
as inferred from numerical simulations, a photometric calibra- 
tion was performed as in previous studies (e.g. Michel-Dansac 
& Wozniak 2004 2006) using the information about the age and 
the metallicity of each particle. I follow closely Michel-Dansac 
& Wozniak (2004). For each stellar particle, given its age and 
metallicity, the mass-to-light ratio T is obtained from a bi-linear 
interpolation into the tables of Bruzual & Chariot (2003), for a 
Salpeter initial mass function (IMF), with mass cut-off at 0. 1 and 
100 M0. Then, luminosity of the particle i (L 1 ), in the V-band, 
is computed according to its mass (M 1 ), birth metallicity (Z^ om ) 
and the time elapsed since its birth (f - fj, om ): 



^v(^ ^born'^born^ 

An age of 5 Gyr is assumed for the initial population at the begin- 
ning of the simulation and a solar metallicity. Other assumptions 
have been tried, in particular imposing an initial metallicity gra- 
dient similar to our Galaxy or setting a younger or older initial 
population, but I found that the results are qualitatively indepen- 
dent of the choice. 

Two-dimensional maps of the age distribution have been 
computed. A low resolution of 41x41 pixels has been chosen 
to smooth out the many meaningless features due to the fluctua- 
tion of the number of particles. For each pixel, the average age is 
the sum of each particle ages weighted either by its mass or by 
its luminosity. These maps are displayed in Fig. Q] and [2] for the 
particular time T = 1.2 Gyr chosen to be illustrative of all others 
snapshots. 

3. Mass- and Luminosity-weighted age distribution 

The ages should be use with caution. Indeed, most likely the 
initial population should have a complex star formation history 
due in part to the assembly of the galactic disk and bulge and in 
other part to its internal dynamical secular evolution. However, 
to limit the number of free parameters, all the initial population 
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Fig. 1. Top panel: large scale V-luminosity-weighted age. The 
field of view is 32 kpc wide. Regions younger than the mean 
are darker and isochronous curves are plotted as dotted lines. 
Regions older than the mean are brighter and isochrones are plot- 
ted as continuous lines. The number of standard deviation from 
the mean labels the isochrones, spaced by 0.25. Bottom panel: 
same as top panel but for a 10 kpc field of view. 



particles have the same age. The amplitude of age differences 
from one place to another in the simulation is thus meaningless 
since it obviously depends on the age of the initial population. 
On the contrary, when a region is dominated by the luminosity 
of a young population, it remains identified as a low-age region 
whatever the age of the initial population because the mass-to- 
light ratio of the young population is much lower than the old 
one. For this reason, I plot in Fig.Q]and[2]the deviation from the 
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Fig. 2. Same as Fig.[TJbut for maw-weighted age. 



mean age as a fraction of the standard deviation. This highlights 
several regions of interest whereas the basic qualitative results 
are thus independent of the assumption made about the old (ini- 
tial) population. 

At large scale (Fig. [TJtop), the two structures with popula- 
tions younger than the mean, are the stellar bar of « 10 kpc 
length and the outer ring (or pseudoring, see discussion in the 
next Section) that encircles it. The youngest regions along the 
ring are shifted by 10 to 20 degrees since star formation is pro- 
voked by the compression in the enrolling gaseous spiral arms 
which are trailing the bar. However, looking at various calibrated 
snapshots between and 3 Gyr, the distribution and morphology 
of the low-age regions in the ring appear very time-dependent. 
Fig- His thus on ly a particular example of such a distribution. 



At the scale of the bar (Fig. [TJ, the young population is 
mainly visible along the major-axis of the bar. Apart from the 
central kpc, that is the site of a regular star formation due to 
the gas inflow, the youngest regions are mainly located between 
roughly 2 and 4 kpc from the centre. On the contrary to the outer 
ring, these low-age regions are persistent, at roughly the same 
place, that is near the end of the stellar bar. The morphology of 
these regions slightly changes with time but not their location. At 
T = 1.2 Gyr (Fig. [TJ, it is noteworthy that their outermost parts 
looks like trailing spiral arms, the rotation being anti-clockwise. 

It is noteworthy that the intermediate region between the bar 
and the ring is depleted from a young and bright population. In 
this region the age of the population is thus quite similar to that 
of the outer disc, beyond the outer ring. The ages deviate by 
roughly +0.5 standard deviation from the map average. 

The luminosity calibration in the V-band obviously enhances 
the region of low-age particles since the mass-to-light ratio of 
these populations is lower than for the initial one (assumed here 
to be 5 Gyr old at the beginning of the simulation, so 6.2 Gyr 
for the snapshot displayed). The maw-weighted is much more 
meaningful for a dynamical analysis (cf. Fig. [2]). 

4. Discussion 

It is noteworthy that the bar low-age regions are just inside or 
close to the ultra-harmonic resonance (or 4/1, hereafter called 
UHR). In Fig. [3] I select all the particles of the young population 
for which 2 < R < 4 kpc and \y\ < 1 kpc at T = 1.2 Gyr. The 
linear UHR radius is ss 4.24 kpc. The formation history of these 
particles is obviously the staking of several episodes of active 
star formation occurring at various places in the galaxy as shown 
by the star formation history of these regions (cf. Fig.[3j. 

The fact that all these particles lie near the UHR along the 
major-axis of the bar does not mean that they were born at this 
place and will remain for a long time. Indeed, the bar dynamics 
is particularly efficient in dissolving these regions since 50 and 
100 Myr later, all these particles are spread out in a wider re- 
gion. Particles trapped by elliptical-like orbits (e.g. x\ family in 
the equatorial plane) remain inside a region approximately lim- 
ited by the UHR radius whereas particles with higher energy can 
fill a wider region limited only by the corotation. For some snap- 
shots, particles have even enough energy to possibly cross over 
the corotation and lie in the disc for a while (the so-called 'hot' 
population). 

The explanation for the persistence of the low-age regions is 
then straightforward. Indeed, particles trapped by elliptical-like 
families of orbits spend most of their time at the apocentre of 
their orbits. The accumulation of young particles on the inner 
side of the UHR is thus simply due to the vanishing radial veloc- 
ity in this region. This is thus a pure effect of the bar dynamics. 

These low-age regions have been already observed in a few 
barred galaxies (Perez et al. 2006; S. Courteau, L. MacArthur 
private communication), using various age and metallicity indi- 
cators. In particular, Perez et al. (f2006) recently claimed that in 
their sample of 6 barred galaxies, the outer parts of all stellar 
bars are younger than the internal parts. To make a crude com- 
parison with their radial plots of Wfi and H3a indices I display in 
Fig.|4]a cut made through the age distribution map of Fig.Q] The 
width of the averaging 'slit' has been taken similar to observa- 
tions, that is 360 pc wide. The age radial profile displays a strong 
gradient as in some observations mentioned above. However, a 
definitive comparison of population ages with Hfi or HSa age 
indicators certainly deserves an in depth analysis using a spec- 
tral calibration of the simulation particles which is beyond the 
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Fig. 3. Particles of the young population belonging to the regions along the stellar bar with the lowest mean age. From left to right: 
distribution of particles selected at 2 < R < 4 kpc and \y\ < 1 kpc at T = 1 .2 Gyr; star formation history displayed as the distribution 
of the epochs of particle formation; particles positions after 50 Myr and 100 Myr. In all spatial plots, the bar has been rotated as to 
be aligned with the x-axis, the red circle is the ultra-harmonic radius, the green ellipse the corotation. 
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Fig. 4. Cut of 360 pc wide along the bar major-axis through the 
V-luminosity-weighted age distribution displayed in Fig. Q] The 
width of the slit is similar to the observations of Perez et al. 2006 
See the text for the explanation of age units. 

scope of this Letter. Moreover, a more detailed comparison with 
these observations will certainly allow some tuning of the initial 
population parameters of my chemodynamical simulations. 

In the central kpc, the situation is quite more complex. 
Following the previous reasoning, the tangential velocity of par- 
ticles on xi-like orbits is maximum at the pericentre. The time 
they spend in the central region is thus very short. However, sev- 
eral generations of particles are created in the central kpc. This 
population is, by the way, responsible of the so-called cr-drop 
phenomenon (Wozniak et al. I2003I I. But, trapped into the po- 
tential well, they mainly remain confined in the circumnuclear 
region (Wozniak & Champavert 2006). As long as the central 
region is forming stars, the mean age will remain several stan- 
dard deviation lower than the average. 

Age gradients and spatial fluctuations are however not ex- 
pected to be very high except inside the bar as discussed above. 
Indeed, for particles not trapped inside the UHR, the gravita- 
tional torques exerted by the bar on the whole disc is a very 
efficient dynamical process to disperse young population. Fig. [3] 
shows that the particles located in the low-age regions at T = 
1.2 Gyr are rapidly scattered. The diffusion timescale is rather 
short, all particles being mixed in less than 100 Myr. This is in 
agreement with the findings of Friedli et al. (119941 1 that the bar 
is responsible to shallow any metallicity gradient outside coro- 
tation. 



The low-age ring that encircles the bar must be undoubtedly 
associated with the so-called outer resonant ring phenomenon, 
since the linear Outer Lindblad Resonance is located at roughly 
8.84 kpc. According to Buta & Combes (1 19961 1, this sort of rings 
appear as blue components, sites of active star formation, on 
deep images, in about one fifth of all spiral disk galaxies. It is 
noteworthy that one third hosts pseudorings that are broken or 
partial rings made up of spiral arms. The feature displayed in 
Fig. [His thus rather a pseudoring than a ring. 

5. Conclusion 

Two diametrically opposed regions, located near the end of stel- 
lar bars are expected to host a composite stellar population with 
a mean age lower than the average. These two low-age regions 
have a dynamical origin since they are due to the accumulation 
of young stellar populations trapped in elliptical-like orbits along 
the bar, near the ultra-harmonic resonance. 

The central region is also a low-age region because of the 
sustained star formation activity generated by the gas accumu- 
lation inside the ILR. Thus, between the nuclear region and the 
ends of a bar, the population age displays strong gradients or 
varying slopes, in agreement with recent observations. 

Another low-age region is the outer ring located near the 
Outer Lindblad Resonance, but the presence, strength and mor- 
phological details of this structure is time-dependent. 
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